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Interplay of simultaneous creation, annihilation, propagation, and relaxation of an excitation in 
molecular condensates interacting with an ultrashort quantum optical pulse is studied in general and 
specialized to a dimer. A microscopic model appropriate for such systems (with a strong exciton- 
phonon coupling) is presented. It also incorporates effects of (quantum) noise in the optical field. 
A variety of new features in the initial stage of excitation dynamics (when it is being created) is 
OS ' revealed; a strong influence of the coherent excitation propagation on the processes of excitation 

On : creation and annihilation in a molecule strongly interacting with phonons is the most remarkable 

one. 
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1. Introduction 

Over the last decades, a great deal of attention has been devoted to the study of dynamics of excitations in systems 
with strong exciton- or electron-phonon interactions like in molecular condensates ( [QJ^J^,^,||,^U^] ; for review, see, 
Refs. [jiP| |l0|JTl| , |l2l] ) . A strong exciton-phonon interaction leads to polaron formation which results in quasicoherent 
propagation of excitations. Quasicoherent propagation has been extensively studied for various models of the exciton- 
phonon interaction based on the (generalized) master equation approach (for review, see, e.g., Ref. U). However, 
all previously developed models assumed that an excitation already exists at the beginning of the interation with 
phonons (e.g., as a result of the interaction with an ultrashort pulse). 

An alternative approach to quasicoherent propagation in molecular condensates is based on the Green functions 
technique |l^,[l4| . But it has in general a different range of validity in comparison with the master equation approach. 
' On the other hand, the optical Bloch equations provide the description of a two-level system with relaxation 
Oh, interacting with a classical optical field Jl5| . This approach is appropriate in many experimental situations (see, e.g., 
> ' Refs. |l^,|l^]). However, it fails in systems interacting strongly with phonons. A generalization of the Bloch equations 
in case of a two-level system interacting strongly with phonons and a classical monochromatic optical field has been 
found in Ref. Jl^] . It has been also shown in Ref. jl8| that field dependent relaxation rates emerge for strong optical 
fields. Also intermolecular coherent transfer influences substantially excitation dynamics in systems interacting with 
optical fields |]l9| , p0| , pl|j22ij23]| . Influence of statistical properties of optical fields in interaction with simple matter 
systems (free and interacting atoms in resonant cavities, atoms in traps, etc.) has been developed in quantum optics 
(see, e.g., Refs. @g||,g|3,|8g]). 

However, the above mentioned theories are not appropriate for the description of excitation dynamics in molecular 
condensates interacting with pulsed optical fields. The paper provides such a theory, i.e. the theory applicable to 
the description of excitation dynamics in systems i) with coherent transfer of excitations, ii) with a strong exciton- 
phonon interaction, and iii) being under the influence of optical fields (e.g., in the form of an ultrashort pulse). It 
thus provides a tool for the investigation of the initial stage of excitation dynamics, i.e. when an excitation is being 
created. Understanding of the excitation dynamics in the initial stage is the main goal of the paper. Moreover, the 
theory encompasses also the influence of (quantum) noise in an optical field. 
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The theory is based on projection operator formalism |50 31 32|j33| , p4| . p5[ . The formalism of the theory is more 
complex in comparison with the previously developed theories, but it is much more general. In fact, it generalizes 
the theories of quasicoherent transfer Q which did not take into account interaction with an optical field and thus 
were not able to describe at least the initial stage of excitation dynamics. However, when the optical field is gone, 
the formalism provides the same excitation dynamics as the older theories. 



* Present address: Joint Laboratory of Optics of Palacky University and Institute of Physics of Academy of Sciences of the 
Czech Republic, 17. listopadu 50, 772 07 Olomouc, Czech Republic 



1 



The general theory is specialized to a dimer interacting with an ultrashort optical pulse. The initial stage of excita- 
tion dynamics is then studied in cases in which processes of excitation creation, transfer, relaxation, and annihilation 
occur at the same timescale, i.e. they mutually compete. In such cases, the previously derived theories cannot be 
applied. Influence of the coherent excitation transfer on the processes of excitation creation and annihilation in a 
molecule strongly interacting with a phonon system is the most important effect. Influence of noise in the optical 
field is also investigated. It is shown, that the theory provides the long-time excitation dynamics in agreement with 
equilibrium statistical physics. 

The theory is primarily developed for the description of excitation dynamics in various molecular condensates. 
However, it can also be applied to other physical systems, e.g., to the dynamics of an excitation in an impurity in a 
crystal. 

Section 2 provides a general theory valid for an arbitrary exciton system. It is specialized to a dimer in Section 3. 
The initial stage of excitation dynamics is investigated in Section 4. Influence of an optical noise is studied in Section 
5. Section 6 is devoted to the long-time behaviour of an excitation. Section 7 summarizes obtained results. Appendix 
contains definitions of coefficients entering equations for the dimer. 

2. Description of a general model 

A system under consideration appropriate for the description of excitation dynamics in molecular condensates 
consists of the exciton (or electron), photon, and phonon subsystems with the exciton-photon and exciton-phonon 
interactions (for details, see Refs. |36|,[37]]). Hamiltonians of the free exciton (H e ), photon (Hf), and phonon (H p h) 
subsystems are given as follows 

He — ^ JmnC^Cri, 

Hf = ^ TiLo K a* K a K , 

K 

H ph =J2 m kbtb k , (1) 

k 

where & m (c m ) means the exciton creation (annihilation) operator at the mth site of the lattice (Frenkel excitons 
are considered), a} K (ax) represents the creation (annihilation) operator of the Kth photon mode and b\. (b^) is 
the creation (annihilation) operator of the fcth phonon mode. The exciton operators and c m obey the Pauli 
commutation relations p8| ; i.e. the operator $ m can be expressed as = d^im, where is the creation operator of 
an electron in the excited state and e m means the annihilation operator of an electron in the ground state at the mth 
site. The photon and phonon operators a) Kl ax, b^, and b^ obey the boson commutation relations. The coefficients 
Jmn describe energies of the free exciton subsystem for m — n and coherent transfer in the exciton subsystem given 
by overlaps of wave functions for m n. The excitonless state |0) is assumed to have zero energy. The symbol lok 
(fife) stands for the frequency of the Kth (kth) mode of the photon (phonon) field. The symbol J^k (52k) means 
summation over all photon (phonon) modes and ^ m denotes summation over all exciton states; h is the reduced 
Planck constant. 

The exciton-photon interaction Hamiltonian H e -f in the rotating wave approximation reads p8| ] 

6 >-f = E huj K F% (a K cl + al K c m ) , (2) 

m,K 



where ujk is a typical photon frequency. The exciton-photon coupling constants F^) are given in the dipole approx- 
imation by 

F$ = -T—^T-TT-— eK ■ (m| P |0> exp (iK ' Ym) ■ (3) 

Tiuj Ko V 2e V ujk m e 

Here e (e < 0) is the charge of electron, m e the mass of electron, p the momentum operator of electron, r rn the mean 
position of the mth molecule in a lattice, |0) describes the excitonless state, |m) the state with one exciton localized 
at the mth site of the lattice and the dot means the scalar product; ek is the polarization vector of the Kth mode of 
the photon field, eo permittivity of vacuum, and V the quantization volume of the electromagnetic field. 
The exciton-phonon interaction is described by the interaction Hamiltonian H e ^ p h in the form |38|,[b| 
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H e _ ph = -L V nn k GTclc m (b k + SO (4) 



N ^ 

m,k 

with only the site-diagonal exciton-phonon coupling included. The dependence of the exciton-phonon coupling 
constants G™ on the site index to and the mode index k is determined according to the type of phonons (optical 
or acoustic) and according to the model of the exciton-phonon interaction; N means the number of phonon modes. 
Hamiltonian (4) describes a deformation of the lattice around a given site m after it was occupied by the exciton 
(polaron effect). 

The excitation dynamics in such a complex system can be conveniently described by the generalized master equations 
for the exciton reduced density matrix p(t). The application of the time-convolutionless formalism ]3l| in connection 
with a time-dependent projector f40|| and the assumption that p{t) lies within the space corresponding to the Hilbert 
space spanned by the excitonless state |0) and by states |m) with one exciton at a given site m ("single-excitation 
approximation" ) result in the following set of equations for the matrix elements of the exciton reduced density matrix 
pit) (for details, see Refs. 

^Pmn(i) — ^ ] JmpPpn{t) + ^ ] JpnPmp{t) ~ ^ ] G mn ^p(t)pp n (t) — ^ ^ G nm p(t)p m p(t) 

p p p p 

-^IrnA^Ppnit) ~ E ^n,p(*)Pmp(*) 
P P 

- iF m (t)p 0n (t) + iF*(t) Pm0 (t) + [/„,„(*) + I*, m (t)] p 00 (t), (5) 



-tX>; (*)/*>(*) + + E KpW M*)> ( 6 ) 

p p S )P 

^po„w = JEvpo p w - X^wm*) - Z^wm*) - E W/W*) 

P P P Z 

+ E [4nW + W*)] ftflW -iD^*(*)ftm(*) +ii ? „*(i)Poo(i)- (7) 
p p 

The matrix element p mm describes the probability that an exciton is in a state localized at the site to and poo 
determines the probability of the system to be in the excitonless state. The nondiagonal elements p mn (to =/= n) 
containing information about the mutual coherence between states localized at sites m and n play the dominant role 
in the description of quasicoherent excitation propagation. The elements po m reflecting a mutual coherence between 
the excitonless state and excited states are important for the description of excitation creation and annihilation. 

The coefficients J mn describe the inner dynamics of the exciton subsystem (coherent transfer) . The time-dependent 
coefficients G mn ^ p (t) and G mtP (t) originating in the exciton-phonon interaction describe effects of polaron formation. 
They renormalizc the coefficients J mn (see Section 6 later). This results in quasicoherent excitation propagation. The 
time-dependent coefficients F m (t) are responsible for excitation creation and annihilation caused by the coherent part 
of the optical field, whereas the coefficients I m , P (t), I m , P (t), and I m ,p(t) reflect effects of the incoherent part of the 
optical field (noise). 

The timc-dcpcndcnt coefficients in Eqs. (5 — 7) are correct to the second power of the exciton-photon or exciton- 
phonon coupling constants. However, this does not mean that the solution of Eqs. (5 — 7) is also restricted to the 
second order. 

The time-dependent coefficients G m „ iP (t) and G m>p (t) are given as follows 

G mn , p (t) = f dT ir E ^ ( G k ~ G k) {M^fc) + 1] CX P (-iflfcr) + 7l B (Ml_fc) exp (iSl- k T)} 
J® k 

xJ2G s - k Hs)(T){p\sr(r), 

S 

G n , p (t) = f ° dr^-Yl ^l G k {M«i*) + 1] exp Hfifcr) + n B (nfl_ k ) exp (iO_ fc r)} 
~' k 

x^G s _ fe (n| S )(r)(p| S )*(r), (8) 
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where the symbol 



n B Wk) = * ; 1 (9) 



denotes the mean value of the number of phonons in the mode K being in the equilibrium state and the correlation 
function (p\s)(t) of the noninteracting exciton subsystem is defined by 

(p\s)(t) = (p\exp(-^H e t\\s). (10) 

The time-dependent coefficients F m (t), I m ,p{t), I m ,p(t), and I m ,p(t) are expressed in the form 

F m (t)=u Ko Pg A(t), 
ImAt) = f dTOj 2 Ka FZ o Y,FL Ko 5N{t,T){p\ S y{t-T), 

I m , P (t) = [ dTu; 2 Ko F^ o J2Fl Ko 5N v (t,r){p\s)*(t-T), 

W*) = f dT^ a FL\ o Y,Fl Ka 6N a (t,T)(p\ s y(t-T), (11) 

Jtn 



where 



Fk = ~z -£Ko ■ Hi>\0) exp (<Ko • r m ) . (12) 



The optical field is assumed to be polarized in the direction €k ■ 
The function 



^>=5Vdb <o * (<)> ' (13) 

describes a classical amplitude of the field. Second moments of the optical field are characterized by the functions 

5N(t,r)= J2 o v ? (Sa K (t)6a* K ,(r))f, 

SN v (t, t) = SN(t, r) + V — % exp [-iw K (t - r)] , 

SN a (t,r) = Yl o y f, = (Sa* K (t)Sa* K ,(r)) f . (14) 

The above introduced correlation functions of the photon field operators are defined as follows: 

(a K {t)) f = Tt f {p f (t )a K (t - t a )} , 
(a K (t)a* K , (r))/ = Try |p/ {t Q )tf K , (r - t )a K (t - t )} , 

(a* K (t)a* K ,(T)) f = Try {p f (t Q )tf K {t - io)a^O - *o)} , 

<5ajc(*) = ««■(*) - <ojf (*)>/• (15) 
A detailed analysis of the above general equations is contained in Refs. |36|J[7]| . 

3. Asymmetric dimer 

The model developed in the previous section is now specialized to a dimer consisting in general of two different 
molecules. We assume that the first (second) molecule has the energy E + 2e (E) and that coherent exciton transfer 
between molecules is described by a real constant J. Hence 
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Jn — E + 2e, J 22 — E, J12 — J21 — J- 

The exciton Hamiltonian (16) can be easily diagonalized by the eigenvectors 

1 



(16) 



1+) 



[J|l) + (- e + A)|2>] 



y/2A(A-e) 
-)= 7=k . =[J|l) + (-£-A)|2)] 



v/2A(AT7) 



(17) 



with the corresponding eigenenergies 



where 



E+ = E + e + A, 
E- =E + e- A, 



(18) 



(19) 



This enables us to calculate the correlation functions (p|s)(i) in Eq. (10) and then to determine the time-dependent 
coefficients entering Eqs. (5 — 7). The resulting equations represent a set of nine coupled differential equations with 
a special structure. They can be conveniently written in the matrix form: 



d 


'fti(t)" 




Ji " 




'fti(t) ' 




' Gi(t) 










n 2 {t)_ 




j 2 _ 




. n 2 (t) _ 




g 2 (t) _ 




n 2 (t) _ 





n 2 {t) 



(20) 



The vectors TZ\(t) and lZ 2 (t) are defined as follows; 



fti(t) 



^2(t) = 



P22(t) 
p r {t) 

Pi{t) 

Poo(t) 

pi r (t) 
Pli(t) 
P2r(t) 
P2i(t) 



(21) 



(22) 



The exciton matrix elements in Eqs. (21) and (22) are expressed as follows: 

1 

2 
1 

2i 



Pr{t) 
Pi(t) 



■ C.C. 



[Pl2(t) - C.C.] 



(23) 



Pjr(t) = ^ <{ />()/(/) oxp 



-i(£? + e)t 



c.c. 



c.c. 



.7 = 1,2, 



(24) 



where c.c. means complex conjugate. 

The matrices J\ and J 2 describe coherent transfer in the free exciton subsystem. Effects of polaron formation as 
a result of the exciton-phonon interaction are contained in the matrices 9i(t) and G 2 (t)- Excitation creation and 
annihilation caused by the coherent part of a photon field is described in the matrices F 2 (t) and T^{t). Finally, the 
matrices Ti(t) and Ti{t) contain the influence of noise in the photon field. Definitions of these matrices as well as the 
description of their role in excitation dynamics are contained in Appendix. We limit ourselves only to the discussion 
of main characteristic features here. 

The matrices on the right-hand side of Eq. (20) have a remarkable block structure. The dynamics of the free 
exciton subsystem is governed by the mutually independent vectors 1Z\ (t) and 1Z 2 (t) as a result of the special type of 
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the exciton Hamiltonian in Eq. (1). Owing to the number of excitation conservation as a consequence of our special 
form of the exciton-phonon coupling (4), the latter interaction changes the dynamics of TZi(t) and 72-2 (t), but it does 
not change their mutual independence. Noise in the photon field acts similarly. The coherent part of an optical field 
introduces mutual coupling between the vectors TZ\{t) and Ti-2{i) and thus leads to an effective exciton generation. 
We note that exciton generation can also be caused by noise in the optical field (for details, see Appendix). 

Interaction with the coherent component of the photon field creates two different kinds of paths leading to excitation 

f k q f1 k 

generation. In the first kind, for molecule 1, an excitation (nonzero pn(t)) emerges along the path poo — ► Pio> Poi ^ 
pn including double interaction of the photon field at this molecule. Terms corresponding to such a path are the 
same as those in the exact equations for a two-level system interacting with a classical deterministic time-dependent 



field 1 37). In this case, the solution of our equations is nonperturbative with respect to the exciton-photon coupling- 



constants. For the second kind of excitation paths, the existence of a surrounding molecule is necessary. An excitation 

F^ F^ j F^ J 

can emerge along the paths p o -*° Pio, Pox Pu,p2i -> pn and p 00 ->° p 2 o, P02 -? P12, P21 -* Pu including double 
interaction with the photon field at molecules 1 and 2 and transfer through J. 

The statistics of an optical field is not limited to a classical noise. Quantum description of the optical field provides 
a tool for investigations of the influence of nonclassical properties of light (e.g. squeezing of vacuum fluctuations). 

Time development of the coefficients originating in exciton-phonon coupling reflects a polaron cloud formation 
around the exciton. The coefficients are practically zero for short times and start to act significantly for longer times 
when they renormalize coefficients in the matrices S\ and J2 ■ Their time development ceases for times comparable with 
the phonon relaxation time tr when they reach asymptotic values. The magnitude of renormalization is proportional 
to the strength of exciton-phonon coupling. Asymptotic values of the renormalization are given by the Debye- Waller 
factor as will be shown in Section 6. That means that our model, although perturbative in exciton-phonon coupling, 
is able to describe correctly effects of a strong exciton-phonon coupling (small polaron formation). 

The above presented equations for a dimer (20) extend Bloch equations in three respects: i) they include transfer 
of an excitation from and to a given molecule, ii) they describe interaction with a phonon reservoir on a "microscopic 
level" (polaron formation), and finally iii) they incorporate effects of (quantum) noise in the optical field. 

4. Initial stage of excitation dynamics 

Based on the above introduced model of the dimer, excitation dynamics is investigated under the conditions when 
times characterizing pulse duration, propagation, and relaxation are comparable. In this case, the model provides 
new interesting results. Firstly, excitation dynamics is studied in a two-level system strongly interacting with phonons 
(Subsection 4.1). In comparison with Ref. [|18j the model is valid also for pulsed light and provides a more accurate 
description of effects stemming from the exciton-phonon interaction. Effects of excitation transfer are studied in a 
dimer in which either one molecule interacts with a pulse (Subsection 4.2) or both molecules interact with a pulse 
simultaneously (Subsection 4.3). Such models are appropriate for a variety of physical situations (e.g. two interacting 
molecules of the same or different kind, a molecule surrounded by an environment with which it can exchange energy, 
etc.). 

Investigation of the excitation dynamics is based on the numerical solution of Eq. (20) for the model of the exciton- 



phonon interaction discussed in Appendix (for details, see, Ref. |37|). The real envelope of the optical field A(t) (for 
definition, see (A9) in Appendix) is assumed in the form A(t) = AA n (t), where A determines the strength of the 
optical field and the function A n {t) has the form 

A n (t) = 1, for t < n 

= exp^-- — — ^ , for t>n. (25) 

Pulse duration is characterized by the constant t\ + t%. The constants F\ = %u)k Fko^ anc ^ ^2 — Twk Fk are 
used in the further discussion. 

We introduce the following simplified notation: Po(t), Pi[t), and P2(t) denote the probabilities that the exciton 
system is in the excitonless state (po(t) = poo(t)), in the state localized at molecule 1 (pi(t) = pn(t)), and in the 
state localized at molecule 2 (p2(t) = P22{t))', p r (t) and Pi(t) mean the real and imaginary parts of pi2(t). Meaning 
of parameters of the system under consideration is schematically shown in Fig. 1. The energy and time scales are 
introduced so that quantities in energy units are in eV and time is in femtoseconds. 



G 







FIG. 1. The scheme of asymmetric dimer; Fi and F2 are the exciton-photon coupling constants at molecules 1 and 2, G\ and 
G2 the exciton-phonon coupling constants at molecules 1 and 2, G is the mean exciton-phonon coupling constant, E + e the 
mean difference of energies between the excitonless state and the excited states, 2e the energy difference between the excited 
states at molecules I and 2, J the transfer integral, u>k the carrying frequency of the optical field and 6' describes the detuning 
of the optical field. The splitting (2e) of excited levels can originate in different kinds of molecules or in the influence of local 
environment of molecules. 



4-1. Two-level system 

A two-level system interacting with an optical field shows Rabi oscillations [[l5| (the stronger the field the higher 
the frequency is). Detuning (5') between the carrying frequency of the pulse and the frequency of the two-level system 
means faster oscillations (with the generalized Rabi frequency) and a lower excitation level (smaller p±(t)). When the 
strength of the optical field is constant the system can return to the excitonless state, but the excitonless state cannot 
be reached in the period of the pulse quenching. 

Interaction with phonons leads to polaron formation around the exciton. This causes a successive diminishing of 
the effective strength of exciton-photon coupling and thus smaller values of the probability p% (t) (it cannot take on 
value 1 for a nonzero exciton-phonon coupling constant G) (see Fig. 2). This diminishing originates in lower effective 
values of the electric dipole moment — e(x) caused by a successive formation of the polaron state from a bare excited 
state and in renormalization of the two- level system energy leading to a greater effective detuning 6' . Strong exciton- 
phonon coupling suppresses deexcitation of the two-level system (see curves B, C, and D in Fig. 2). Effects of polaron 
formation are more pronounced for higher temperatures, i.e. for greater mean numbers of equilibrium phonons «b- 




200 400 600 



FIG. 2. The influence of the exciton-phonon coupling constant G on the probability pi (i); G = [A], G = 0.004 [B], G = 0.01 
[C], G = 0.02 [D]; J = 10~ 8 , W = 0, Fi = 0.01, F 2 = 0, e = 0, n = 100, r 2 = 100, Gi/G = 1 + 0.25i, G 2 /G = 1 - 0.25i, 

n B = 0, nn ph = 0.01, hy ph = o.ooi. 



x The magnitude of exciton-photon coupling constants in energy units can be estimated from the mean value of the exciton- 
photon interaction Hamiltonian H e -f as follows: (ff e _/) ~ ^-(p)A ~ j- ([H e , x]) ^ w i(£<= ~ ^a)(^) ^ ~ f(^ e ~ ^s) 123 !^' 
where A is the vector potential of a classical field, E the intensity of the field, u>o the mean frequency of the field, p (x) means 
the momentum (position) operator of electron, £ e (£ 9 ) denotes the energy of the excited (ground) state, and as is the Bohr 

radius; E = 10 7 Vrn -1 , A = 600 nm, and £ e - £ g = 2 eV implies (# e _/) w 5.10" 4 eV. 
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4-2. Dimer with one molecule interacting with a pulse 



Only molecule 1 is assumed to interact with the pulse as, e.g., a consequence of a special orientation of the transition 
dipolc momentum of molecule 2. In case of exciton dynamics, the "single-excitation" approximation requires that 
values of the whole excitation probability of the dimer pi(t) + P2(t) are approximately less than 0.6 (for details, see 
Subsection 4.3). The obtained results are appropriate also for electron transfer as a result of the "single-excitation" 
approximation which rules out differencies originating in Fermi and Pauli commutation relations. In this case, the 
above limitation does not apply. 

4-2.1. Energetically balanced dimer 

The effect of coherent transfer (described by J) on the dynamics of the dimer noninteracting with phonons is shown 
in Figs. 3a, b. The increase of J leads to a faster exchange of excitation between molecules 1 (pi(t)) and 2 (p2{t)) 
(its frequency is given by J). The transfer strongly affects the processes of excitation creation and annihilation at 
molecule 1. Fig. 3a shows that the value of J affects the time in which the excitation of molecule 1 is replaced by its 
deexcitation; especially greater values of J mean an earlier time of deexcitation. This results in smaller values of the 
whole excitation probability pi(t) +Pi{t) for greater values of J. 



0.8 n 0.5 -i 




FIG. 3. Increasing values of the transfer integral J influence values of the probabilities pi(t) (a) and pi(t) (b); J = 10 8 
[A], J = 0.0005 [B], J — 0.001 [C], J = 0.002 [D]; Fi = 0.0005, F 2 = 0, s = 0, US' = 0, n = 1000, r 2 = 200, G = 0. 

Molecule 2 can be excited (p2 ^ 0) even in the time when molecule 1 is not excited (pi = 0) (see Fig. 4 for t w 300). 
In this time molecule 2 is completely decoupled because p r = pi = 0. Such an effect has already been observed in 
Ref. The exciton subsystem can also return to the excitonless state (po = 1, see Fig. 4 for t w 600). However, 
when the pulse becomes weak the exciton subsystem cannot reach the state with po = 1 (see Figs. 3a, b). 




-o.oi 

500 1000 1500 2000 

FIG. 4. A typical time development of the probabilities pi(t) [A] and p^{t) [B] and the nondiagonal matrix elements p r {t) 
[C] and pi(t) [D] in the phononless dimer; Ft = 0.0005, F 2 = 0, J = 0.007, e = 0, M' = 0, t\ = 1000, r 2 = 200, G = 0. 

Interaction with phonons influences the dynamics as follows. In case of small J a decrease of the whole excitation 
probability pi(t) + P2{t) (see Figs. 5a, b) with the increase of G occurs as a consequence of the prevailing effect of 
polaron formation at molecule 1. An increase of the whole excitation probability pi(t) + P2(t) with the increase of 
G for greater J is observed, because J is renormalized by the interaction to its smaller values which suppresses the 
destructive effect of J on the excitation creation. But when G is great enough (according to the value of J) the effect 
of polaron formation at molecule 1 decreases the whole excitation probability. The magnitude of G also influences 
the time when the excitation creation is replaced by its annihilation. Especially, annihilation occurs earlier when the 
polaron formation at molecule 1 prevails for smaller J (compare curves A, B, and C in Figs. 5a, b). 
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FIG. 5. Increasing values of the exciton-phonon coupling constant G mean the decrease of the probabilities pi(t) (a) and 
p 2 (t) (b); G = [A], G = 0.003 [B], G = 0.005 [C]; F 1 = 0.0005, F 2 = 0, J = 0.0005, e = 0, M' = 0, n = 1000, r 2 = 200, 
Gi/G = 1 + 0.25i, G 2 /G = 1 - 0.25i, Tifiph = 0.01, Tvv, = 0.001, n B = 0. 

After the pulse is gone, the frequency of excitation exchange between molecules 1 and 2 decreases with increasing 
G. This manifests the renormalization of J. Also a complete deexcitation of molecules {p\ = p 2 = 0) cannot be 
reached owing to the polaron effect. 

4-2.2. Energetically unbalanced dimer 

The energy difference 2e between the excited states of molecules 1 and 2 affects the phononless system as follows. 
Greater values of the probabilities pi(t) (see Fig. 6a) and pi(t) + P2(t) for greater values of |e| have their origin in 
decoupling of molecule 2 from molecule 1 which partially suppresses the effect of coherent transfer (J). Decoupling 
of molecule 2 then causes smaller values of P2{t) (see Fig. 6b). The case when molecule 2 has a lower energy than 
molecule 1 (e > 0, 6' < 0) cannot be distinguished from the case in which molecule 2 has a higher energy than 
molecule 1 (e < 0, 8' > 0); time development of Po(t), Pi(t), and P2(t) is the same. 




FIG. 6. The increase of the energy difference 2e suppresses the effect of J, this is demonstrated in the time development 
of the probabilities Pl (t) (a) and p 2 {t) (b); e = 0, hS' = [A], s = 0.0005, W = -0.0005 [B], e = 0.001, W = -0.001 [C], 
e = 0.002, W = -0.002 [D]; Fi = 0.0005, F 2 = 0, J = 0.002, n = 1000, r 2 = 200, G = 0. 

The increase of e may result in the decrease of the probability pi (t) +p 2 (t) when exciton-phonon coupling is nonzero. 
Also the above mentioned indistinguishability is lost. 

If molecule 2 has a lower energy than molecule 1 (positive e), values of the probabilities Pi(t) (see Fig. 7a) and 
Pi(t) +P2{t) (compare Figs. 7a,b) increase with the increase of G. This indicates the increasing renormalization of J. 
However, when values of G are too great, then polaron formation at molecule 1 enters into play and the lowering of 
values of pi(t) and pi{t) + P2{t) is observed. 

Increasing positive s leads to greater values of p\ (t) + p 2 (t) in cases when the polaron formation at molecule 1 
does not play an important role (for smaller values of G). This is caused by the decrease of the destructive effect 
of transfer on the excitation creation. However, the increase of positive e means the decrease of p\ it) + p 2 (t) if the 
polaron formation at molecule 1 plays a dominant role. In this case, the destructive effect of transfer on the polaron 
formation at molecule 1 becomes weaker with increasing e. 
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FIG. 7. Increasing values of the exciton-phonon coupling constant G suppress the effect of J and that means the increase of 
the probabilities P i(t) (a) and p 2 (t) (b); G = [A], G = 0.003 [B], G = 0.005 [C]; F ± = 0.0005, F 2 = 0, J = 0.005, e = 0.004, 
M' = -0.004, n = 1000, t 2 = 200, Gi/G = 1 + 0.25i, G 2 /G = 1 - 0.25i, hn ph = 0.01, hj ph = 0.001, n B = 0. 

If molecule 2 has a higher energy than molecule 1 (negative e), the increase of G results in the decrease of the 
probabilities Pi(t) and pi(t) + P2(t) in all cases. Thus the interaction with phonons supports the destructive effect of 
J on the excitation creation. 

The increase of J need not mean only smaller values of pi (t) caused by the destructive effect of J on the excitation 
creation. Especially, when a strong interaction with phonons forms a polaron at molecule 1 suppressing its excitation, 
greater values of J can partially break the polaron formation and admit greater values of p\ (i) (this effect is well 
pronounced for greater positive e). 

An exciton moving on the dimer can be "partially localized" || by the interaction with phonons (see Fig. 8). A 
complete localization of the exciton cannot be reached because coherent transfer cannot be completely suppressed 
(see asymptotic expressions in Section 6). 




FIG. 8. The "partial localization" of exciton in the dimer, described by the probabilities pi(t) [A] and p2(t) [B] and by the 
nondiagonal matrix elements p r (t) [C] and p t (t) [D]; Fi = 0.0005, F 2 = 0, J = 0.0005, e = 0, HS' = 0, n = 1000, r 2 = 200, 
G = 0.005, Gi/G = 1 + 0.25i, G 2 /G = 1 - 0.25i, m ph = 0.01, h lph = 0.01, n B = 0. 

4-3. Dimer with both molecules interacting with a pulse 

We first address the validity of "single-excitation" approximation considering a two-level system and a symmetric 
dimer with J being practically zero and comparing their levels of excitation. Omission of the two-exciton state in 
the dimer manifests itself in lower values of the excitation probability p\{t) (or P2(t)) in comparison with that for a 
two-level system. The deviations can be estimated for a given value of the excitation probability from curves in Fig. 9. 
In general, the "single-excitation" approximation is very well accepted for values of the whole excitation probability 
Pi{t) + P2(t) up to w 0.2. The deviations in pi(t) + p^(t) are less than about 10 % for values of pi(t) + p-z(t) up to 
w 0.6, which is also well acceptable. 

Results of this section are not applicable to the electron dynamics because only one excitationless state for the 
whole system has been considered. 
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FIG. 9. The time development of the probability p\(t) shows the accuracy of the "single-excitation" approximation; F 2 = 
[A], F 2 = 0.0002 [B]; Fi = 0.0002, J = 10~ 8 , e = 0, W = 0, n = 1000, r 2 = 200, G = 0. 

4-3.1. Energetically balanced dimer 

Coherent transfer (J) affects the excitation dynamics similarly as it is discussed in Subsection 4.2; i.e. the increase 
of J in the phononless system means the decrease of the whole excitation probability p\{t) + P2(t). The exciton 
subsystem can return to the excitonless state (po = 1). 

The increasing interaction with phonons (G) decreases pi(t) +P2{t) for small J as a result of the polaron formation 
at both molecules. However, the increase of G leads to the increase of Pi(t) + P2(t) for greater values of J because 
the destructive effect of transfer on the excitation creation is partially suppressed. 

4-3.2. Energetically unbalanced dimer 

We assume that molecule 1 is pumped resonantly and molecule 2 nonresonantly in the further discussion. Despite 
the effect of coherent transfer, small values of e lead to a strong asymmetry in the time development of the probabilities 
Pi(t) and P2(t) (see Fig. 10) having its origin in the nonresonant interaction of molecule 2 with the pulse. 

0.08 




r>C0 



1000 

t 



1500 



2000 



FIG. 10. A strong asymmetry in the time development of the probabilities pi(t) [A] and P2(t) [B] is caused by nonzero e; 
Fi = 0.0002, F 2 = 0.0002, J = 0.002, e = 0.0005, W = -0.0005, n = 1000, r 2 = 200, G = 0. 

The energy position of molecule 2 with respect to molecule 1 (the sign of e) strongly influences the excitation 
dynamics (see Fig. 11); the probability p\ (t) reaches greater values for molecule 2 having a lower energy than molecule 
1. Increasing |e| smooths this asymmetry. 



0.00 




FIG. 11. The sign of the energy difference 2e strongly influences the time development of the probability pi(t); e = 0.002, 
W = -0.002 [A], e = -0.002, W = 0.002 [B]; Fi = 0.0002, F 2 = 0.0002, J = 0.002, n = 1000, r 2 = 200, G = 0. 
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Interaction with phonons modifies the excitation dynamics as follows. For molecule 2 having a lower energy than 
molecule 1 (e > 0), the probability pi(t) + P2(t) increases with increasing G until the range of values of G is reached 
in which the polaron formation at molecule 1 causes the decrease of pi (t) + p2 (i) . The higher the e is the less the 
destructive effect of J on the polaron formation at molecule 1 is. The excitation dynamics is then influenced for 
weaker G. 

When the energy of molecule 2 is slightly higher than that of molecule 1 (e < 0), the increase of G leads to the 
decrease of pi(t) (effects of transfer and polaron formation support each other). But an increase of P2(t) occurs 
because interaction with phonons tunes molecule 2 to the resonance with the field. The probability p\{t) + P2(t) also 
increases with increasing G. For greater values of |e|, the probability pi(t) +P2(t) decreases with increasing G. 

After the pulse is gone the excitation dynamics is the same as that discussed in Subsection 4.2. 



The obtained numerical results elucidate the validity of perturbation approximation. We consider the perturba- 
tion approximation in the exciton-photon coupling constants to be reliable for arbitrary times at least for classical 
deterministic fields. The positive semidefiniteness of the exciton reduced density matrix has been slightly breaked for 
longer times and greater values of G. This clearly indicates that it is caused by the omission of terms proportional to 
the third and higher powers of the exciton-phonon coupling constants. However, deviations are small for real values 
of parameters and the perturbation approximation in the exciton-phonon coupling constants is also justified. 

5. Influence of optical field fluctuations 

Effects originating both in the amplitude and phase fluctuations of a pulsed optical field on the excitation dynamics 
are studied. 

We assume that the noisy part of the photon field amplitude is proportional to the strength of its coherent part. 
The normalized photon field correlation function 5N n (t,T) has then the form (the strength of the field is absorbed 
into Fi and F2): 

SN n (t, T ) = A n (t)A n (T)(5A n (t)6A* n {T)) f , (26) 

where the normalized envelope of the pulse is given in Eq. (25). The optical field fluctuations are described by the 
amplitude deviation 5A n (t), 



normalized with respect to the envelope A(t) defined in (A9) in Appendix. A stochastic model ( |26|, p. 137) for both 
the amplitude and phase fluctuations of SA n (t) provides the relation 

(SA n (t)SA n (r)) f = (\6An\ 2 ) f exp[-i(uj Ko +LJ s )(t-T)] (exp [icf>{t) - #(r)])/. (28) 

The moment (|<L4 ra | 2 )/ involves averaging over amplitude fluctuations, uj s denotes the frequency shift of amplitude 
fluctuations and the factor (exp [i(f>(t) — i</>(r)]) / involves averaging over the phase <f>{t). E <j>(t) - 4>(t) = J'gIt' Aw(t'), 
where Alo(t) represents the Gaussian Markovian process, then Eq. (28) can be rewritten into the form 

(SAntySA^fj)) f = n s exp [-i(u Ko + u s )(t - t)} exp [- 7a |t - r|] . (29) 

The quantity n s = (|<5.4 n | 2 ) / characterizes the strength of noise, whereas j s describes the strength of phase correlations 
«Aw(i)Aw(i')) = 2 ls 5(t -t')). 

The correlation function 8N n (t,T) in (26) then gets the form: 

5N n (t,r) =A\ l (t)Al(T)n s eiq>[-i(uj Ko +uj s )(t-T)]e-x.p[-j s \t-T\]. (30) 

We further omit vacuum fluctuations in the photon field and thus SN v (t,r) = 5N(t,r) (see Eq. (14)). We also 
assume SN a (t,r) = as a consequence of averaging over the phases of 5a* K {t) and 5a* K , (r) in (Sa* K (t)5a* K , (r))f in 
Eq. (14). 

The influence of noise on the excitation dynamics is described by the matrices T\{t) and T^t) given in Eq. (A7) 
in Appendix. The time-dependent coefficients Mi(t) = Mi(t), . . . , = P2(t) defined in Eqs. (A10 — A12) contain 
the normalized correlation function <5iV„(i,r) given in Eq. (30) instead of 6N(t,r). 

Optical noise influences the dynamics of a two-level system in such a way that it forces the system to evolve in the 
direction to the state with po = p\ = 1/2, pio = (see Fig. 12). The increase of n s leads to the increased influence 
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of noise. The decrease of j s means the increase of the correlation time of noise t s (t s — l/j s ) and results in the 
increased effect of noise. 
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FIG. 12. A typical excitation dynamics in the two-level system interacting with a noisy pulse described by the probabilities 
po(t) [A] and pi(t) [B]; Fi = 0.01, F 2 = 0, J = 10~ s , e = 0, W = 0, n = 1000, r 2 = 200, G = 0, n s = 0.1, 7 S = 0.01, lo s = 0. 

For values of the exciton- photon coupling constants F\ and F 2 and pulse durations (n + r 2 ) used in Section 4, the 
influence of noise with reasonable values of parameters (n s < 0.1, 7 S > 0.01) is negligible. 

These results show that noise in ultrashort pulsed optical fields does not influence substantially experimental results 
under standard conditions. Effects originating in nonclassical properties of photon fields are the matter of further 
investigations. 

6. Long-time behaviour 

Eq. (20) is analyzed for times when the photon field (the pulse) does not act on the exciton subsystem and for 
times longer than the relaxation time tr of the phonon reservoir. Vacuum fluctuations of the photon field are also 
omitted and hence, e.g., the exciton decay processes are not taken into account. Then the third matrix in Eq. (20) is 
zero and the sets of equations for 1Z\(t) and IZ2 (t) are independent. The time-dependent coefficients in the matrices 
Gi(t) and G2 (t) can be replaced by their asymptotic values for t — to — ► oo. 

The excitation dynamics is then driven by the following equations: 





' Pu{t)' 




d 


P22(*) 


1 


dt 


Pr(t) 


~ h 




I Pi(t) . 





-2J 
2J 
2e 
J -J -2e 



~Pn(t)- 




- 








P22(t) 













Pr(t) 






Cas 


Eas 


I Pi(t) . 




. -Bas 


-D as 








-F a ., 





~Pu(t)- 




P22{t) 




Pr(t) 




L pi(t) J 



(31) 



The coefficients A as , . . . , F as express asymptotic values of the time-dependent coefficients A(t), . . . , F(t) given in (A3) 
in Appendix for t — t — > oo; 



ttJ 



%N 
J 

fiN 

ttJ 
hN 



\G\ - G 2 k \ 2 {e [2n B {hn k ) + 1] + A} 5{hfl k - 2A), 



G\\ 2 V 



n 2 k 

Q 2 k - 4A' 2 



^4as 

C as = ^J2 l G fe - G fe! 2 [ 2 "B(^ fe ) + 1] + A} 5(m k - 2A), 
D as 

E as 
F 

± as 



n 2 - 4A' 2 



k 

2ttJ 2 
hN 



2n B (^ fc) + l] + J|}, 



[2nB (Ml fc ) + l]-^} 



J2 \G\ - G\\ 2 [2n B (hn k ) + 1] S(hn k - 2A), 



h E ( G fc ~ G fe) ( G -fe + G -fe) fife ' 



(32) 



T 3 ' denotes principal value. 

For s = (it implies that £> as = D as ) and J <C AW p h (AW p h being the width of phonon spectrum) the transfer 
integral J at the positions (4,1) and (4,2) in Eq. (31) is renormalized to a new value 



13 



J ren = J-hB as = Jexp(-2VK) + o(G 2 ), 
1 
27V 



W -.J-Y^lGl-GlfpriBm^ + l}, (33) 



where W denotes the Debye- Waller factor describing the small polaron formation p^j44[ |. 

The transfer integral J at the positions (4,1) and (4,2) is renormalized twice and J at the positions (1,4) and (2,4) 
is not renormalized. But the effect of renormalization on the matrix elements pn(t) and P22(t) is correct as can be 
shown from the equations for pu(t) and p22(t) with the nondiagonal matrix elements puft) and p2i(t) excluded [Q. 

The asymptotic state of the dimer (i.e. the state reached for t — to — > oo) has the form correct up to G 2 as follows, 

(pTi,P22,Pr,pT)^ (-7iA + 72£,-7iA - 72 s,J l2 ,0), (34) 

where 

j- \g\ - G\\ 2 [2n B (hn k ) + 1] 5(nn k - 2A), 

k 

^Ei G fe- G "i 2(5 ( M ^- 2A )- ( 35 ) 



71 N 

k 



72 N 

k 



Its form written in the basis of the eigenstates |+) and |— ) (17) 

_ A - e 



pas _ A + £ 



exp(-2/3A), P 5f_ = p™, = 0, (36) 



is in agreement with the result of equilibrium statistical physics for an already created excitation. This indicates that 
the theory is valid also for longer times despite the perturbation approximation applied in the derivation. 

7. Conclusions 

A microscopic theory of the excitation dynamics in molecular condensates interacting with optical fields (with 
ultrashort pulses) has been developed. It has provided a tool for the investigation of the initial stage of excitation 
dynamics, i.e. when an excitation is being created. 

The theory has been applied to a dimer under the conditions when times characterizing pulse duration and propaga- 
tion and relaxation of an excitation are comparable. The following results have been obtained. In this case, coherent 
transfer strongly influences the processes of excitation creation and annihilation. The stronger the coherent transfer 
is, the lower the excitation level is. The influence of coherent transfer is suppressed by the increase of energy difference 
between excited states. Interaction with phonons leads in general to polaron formation. This leads to lower excitation 
levels on one side, but on the other side it partially preserves the excitation from being annihilated. Interesting 
behaviour occurs when effects of coherent transfer and polaron formation compete. Interaction with phonons can 
partially suppress the destructive effect of coherent transfer on the excitation creation; this leads to higher excitation 
levels. This occurs when the energy of a molecule being excited is higher than that of surrounding molecules. On the 
other hand, interaction with phonons can support the effect of coherent transfer on the excitation creation and thus 
can lead to the further lowering of the excitation level. This is the case when the energy of a molecule being excited 
is lower than that of surrounding molecules. The theory also provides a long-time dynamics being in agreement with 
equilibrium statistical physics. 

It has been shown that noise in an ultrashort pulse does not influence substantially the initial stage of excitation 
dynamics under standard conditions. However, the theory is also suitable for the description of effects originating in 
nonclassical properties of interacting optical fields (e.g., those with squeezed vacuum fluctuations). 
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APPENDIX A 



Definitions of matrices in Eq. (20) 

Appendix contains the definitions of the matrices Ji, J2, Gi(t), Q2(t), T\{t), ^(t), ^(t), and !Fi(t) in Eq. (20). 
The influence of their elements on the excitation dynamics is also discussed. 
The matrices J7i and J2 describe the dynamics of the free exciton system: 



-2J 
2J 
2e 
J -J -2e 




n 



-e -J 
e J 
-J e 
J -e 



(Al) 



The matrices Q\{t) and G2(t) stemm from the exciton-phonon interaction: 



Gi(t) 







A C E —F 

B -D F E 





&(*) 



A l —B\ & -£>i 

Si Ax d x a 

_ i?2 C2 -B2 A2 



(A2) 



The time-dependent coefflcients A(t), . . . ,D 2 (t) are given as follows 
A(t) 



The functions 



C(t 
D(t 

E(t 
F(t 
A^t 

Bi(t 

A 2 (t 

B 2 (t 

Ci(t 

Dx{t 

C 2 (t 

D 2 (t 



Je 
A 2 



-^G 2 .92, 2 -^G 2 53,i, 

t/ £ ^2 — ^ /t2 — 

72 G .92,1 + ^T-G 33,2, 



A 2 
Je 
A 2 



2A 
2A' 



2- j2 2- 

G 51,1 — 32,1, 

-G 2 3i, 3 , 

J 2 

G 2 gi,4 + ^G 2 [-g 2 ,4 + 52,5 + 52,10] , 

J 2 

G 2 3i,s + ^G 2 [-52,6 - ff2,8 + 32,9] , 

J 2 

G 2 #i,7 + ^G 2 [.92,5 - 32,7 - 52,10] , 

J 2 

G 2 5i,n + ^G 2 [52,6 + 52,9 - 52,11] , 



Je 
A 2 



G 2 [52,4 ~ 52,5 - 52,10] 



J 

2A 



G 2 [53,6 + 53,8 - 53,! 



Je J 

^■G 2 [52,6 + 52,8 ~ 52,9] + ^-G 2 [-5 3j4 + 53,5 + 53 

^■G 2 [52,5 - 52,7 - 52,10] + ^-G 2 [-53 



10J 



Je J 
-^G 2 [52,6 + 52,9 - 52,11] + [53, 



6 — 53,9 + 53,ll] 
5 — 53,7 — 53,lo] • 



9i,j{t) 



t-t 



(A3) 



(A4) 
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t-t 



g 2J (t,A')= / dr^r) sin 2 (A'r) 



53,i(t,A') 



drgj^r) sin(2A'r), 



J = 1, • 
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(A5) 



describe the response of the exciton subsystem to the phonon one. The phonon subsystem is characterized by the 
functions □ 



91 (t 

32 (t 

33 (r 

34 (t 

35 (t 

36 (t 

37 (t 

38 (t 

39 (t 
3io (t 
511 (t 



i- £ ^-|G^ - G|j 2 [2» B (ftfifc) + 1] cos(fi fc r), 



G 2 N 



G 2 N 



\- Y, VfoGl - G 2 k ){G l _ k + G 2 _ k ) sin(^r), 



G 2 N 



\- ]T Vi\Gl\ 2 [2n B (hn k ) + 1] cos(fi fc r), 



G 2 N 



i- Y ^lMG\G 2 _ k ] [2n B (hn k ) + 1] cos(tt fc r) 



G 2 N 



\- J2 n 2 lm[GlG 2 _ k ] [2n B (hn k ) + 1] cos^r), 



G 2 N 



l — Y n 2 k \G 2 \ 2 [2n B (hn k ) + 1] cos(fi fc r), 



G 2 N 



\ 1 Y^t\Gl\ 2 sMQ k T), 



G 2 N 



l-Y^lMGlGl^i^r), 



G 2 N 



l~Y^iMGlG 2 _ k ] S m(^kT), 



G 2 N 



L-YnHGH 2 S m(Q k T) 



G 2 N 



(A6) 



The constant G having the meaning of the mean exciton-phonon interaction constant has been introduced into 
Eqs. (A3) and (A4) as well as into the definitions given in Eq. (A6) in order to get gi(r), . . . , gn(r) dependent only 
on the dispersion of coupling constants. The symbols Re and Im denote real and imaginary parts. The new symbol 
A' = A/h has been introduced here. 

The time-dependent coefficients B(t) and D{t) given in Eq. (A3) renormalize the transfer integral J. The coefficients 
A(t) and C(t) are important for relaxation to equilibrium state (see, e.g., Refs. El 42 4^]). 

The matrices !Fi(t), T 2 {t), ^(t), and F±(t) originate in the interaction with an optical field: 



' 2M 1 





20i 


20 2 


-2Mi - 





2Nx 


2Pi 


-2P 2 


-2Ni 


Pi 


Oi 


Mi + Nx 


-M 2 + N 2 


-Ox - Pi 


-Pi 


2 


M 2 - N 2 


Mi + Nx 


p 2 -o 2 


. -2Mi 


-2JVi 


-2Px - 20x 


2P 2 - 2Q 2 


2Mx + 2Ni . 



2 In numerical calculations, we assume that hQ, k G k = Gi in one half of the fc-space and hQ k G k = G* in the remaining half of 
the fc-space (i = 1, 2) (for details, see, Ref. J37|). This assumption is in agreement with hermiticity of H c - P h- Further, the mean 
numbers of phonons n-B{Kl k ) axe assumed to be fc-independent (nB{hQ ko ) — Ub)- The remaining summations ^ fc sin(fifcT) 
and cos(Q k r) in Eqs. (A6) are replaced by the expressions sin(fl p hT) exp( — 7 P hr) and cos(f2 p fcT) exp(— 7 P hr), respectively. 

The frequency Q. p h then characterizes a mean phonon oscillation frequency and 7 P h describes damping originating in dephasing. 



1G 



W) = 



-F 4 (f) 



2K X -2K 2 
2Li 

— L2 — il ^2 





Kx 



-2ifi 2^2 -2Li 2L 2 

— isTi -Li L 2 #1 
^2 L 2 Li -if 2 
-L x -Kx -K 2 Lx 
L 2 K 2 -Kx -L 2 

2Mj + iVx - 2Mx 2M 2 + N 2 - 2M 2 

-2M 2 -N 2 - 2M 2 2Mx + Nx + 2M X 
Px-Ox- Px P 2 -0 2 - P 2 
-P2 -0 2 -P 2 Px+Ox+ Px 



Ox -Ox- Px 2 -0 2 - P 2 

-0 2 -6 2 -P 2 Ox+Ox+ Px 

2Nx + Mx - 2Nx 2N 2 + M 2 - 2N 2 

-2N 2 - M 2 - 2N 2 2Nx + Mx + 2N ± 



(A7) 



The coefficients Kx(t), K 2 (t), Lx(t), and L 2 {t) describing the influence of the coherent part of an optical field have 
the form (the constants and are assumed to be real): 

Kx(t) = -uj Ko F^Iux [A{t)exp(i6't) 
K 2 (t) = uj Ka Fk Re [A(t)exp(i5't) 
Lx(t) = -oj Kq FI q Ixxx [^(i)exp(^'t) j 

L 2 {t) = u Ka F 2 Ko Rc \A{t) exp(i<y'i)] . (A8) 

The symbol 5' denotes the frequency mismatch (6' = (E + e)/% — u>k ) and the envelope A(t) of the field is defined 
as follows, 



A(t) = A(t) exp(ioj Ko t). 



(A9) 



The coefficients Mx(t), . . . ,P 2 (t) reflect statistical properties of the optical field (noise) and can be expressed in the 
form 



Mx(t) = (F^fix - (Fk ) 2 ^U - F^F^u, 




M 2 (t) = (Fk )\ + (Fk ) 2 ^ 2 + Fk F 2 Ko ^ 2 , 




Nx(t) = (F 2 Ko )^x + (F 2 Ko f^u - Fk P£ ±U, 




N 2 (t) = (F 2 K J 2 i 3 -(F 2 Ko ) 2 ^i 2 + F^F 2 Ko ^i 2 , 




Ox(t) = -(Fk f^U + Fk F 2 Ko ix + FkA^u, 




2 (t) = (Fk ) 2 ^i 2 + F' Ko F 2 Ko i 3 -F' Ko F 2 Ko ^i 2 , 




W) = -(F 2 Ka ?^x + Fk F 2 Ko ix - FkA^i 4 , 




P*{t) = (Flf^ 2 + Fk F 2 Kol3 + Fk F 2 Ko ^ 2 . 


(A10) 


The functions ix(t, A', 5'), . . . , ii{t, A', 5') characterize the response of the exciton 


subsystem to the photon field: 



ix(t,A',5')=w 2 Ko [ dr cos[A'(i- r)] Re 5N(t, r) exp [iS'(t - r)] 

J t 

i 2 (t, A', 5') = lo 2 Ko [ dr sin [A'(t - r)] Re \SN(t, r) exp [i5'(t - r)] 
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i 3 (t, A', S') = uj 2 Kq f dr cos [A'(t - t)] Im \dN(t, r) cxp [iS'(t - r)] 

i 4 (t, A', <5') = / dr sin [A'(t - r)] Im UiV(t, r) exp [iS'(t - r)] 

The photon field correlation function 5N(t, r) is of the form: 

5N(t, t) = 5N(t, t) exp [iu Ko (t - r)] . 



(All) 



(A12) 



The coefficients Mi(i), . . . , P2W are defined similarly as the coefficients Mi(i), . . . , Pz^) m Eqs. (A10) and (All); 
only the photon field correlation function 



5N v (t, t) = SN v (t, t) exp [iu Ko {t - t)} 



(A13) 



occurs in Eq. (All) instead of SN(t, r). Thus, the coefficients with bars include in addition effects of vacuum 
fluctuations. 

Also the coefficients Mi(t), . . . , P2(t) are defined similarly as the coefficients M\(t), . . . , P2(t) in Eqs. (A10) 
and (All): only the expression SN(t, r) exp [i5'(t — r)] in Eq. (All) must be replaced by the expression 
SN a (t 7 t) exp [— i5'(t + r)], where 



SN a (t, t) = SN a (t, t) exp [-iuj Ko (t + r)] 



(A14) 



The coefficients 02(t) and P2(t) rcnormalize the transfer integral J both at the positions (4,1), (1,4) and (4,2), 
(2,4) in the matrix J7i, in contrast to the coefficients originating in the exciton-phonon interaction. 

The influence of the coefficient M\{t) (N\(t)) in the equations for p\ r {t) and pn(t) (p2r(t) and P2i(t)) is remarkable. 
When, e.g., M\(t) (N\(t)) is negative, it represents damping of p\ r (t) (p2r(t)), but at the same time amplification of 
Pu{t) (p2i(t))- This property is connected with phase relations in the photon field reflected by SN a (t,T) given in Eq. 
(A14). 
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